Steady-state conduction in self-similar billiards 
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The self-similar Lorentz billiard channel is a spatially extended deterministic dynamical system 
which consists of an infinite one-dimensional sequence of cells whose sizes increase monotonously 
according to their indices. This special geometry induces a non-equilibrium stationary state with 
particles flowing steadily from the small to the large scales. The corresponding invariant measure 
has fractal properties reflected by the phase-space contraction rate of the dynamics restricted to a 
single cell with appropriate boundary conditions. In the near-equilibrium limit, we find numerical 
agreement between this quantity and the entropy production rate as specified by thermodynamics. 
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Over the last two decades, the non-equilibrium sta- 
tistical mechanics of chaotic dynamical systems has re- 
ceived a great deal of attention, in particular regarding 
the positivity of entropy production and its connection 
to dynamical properties characteristic of chaos, such as 
the Lyapunov exponents P, 0, H, Hi ■ 

Of particular interest is the two-dimensional periodic 
Lorentz gas with external forcing and Gaussian iso- 
kinetic thermostat (GIKLG) This is an elementary 
model of electronic conduction under a non-equilibrium 
constraint. The thermostat is actually a mechanical con- 
straint and acts so as to remove the energy input from the 
external forcing 0]. Under the action of this thermostat, 
the kinetic energy remains constant and thus fixes the 
temperature of the system; no interaction with a hypo- 
thetical environment is needed in order to achieve ther- 
malization. Rather, dissipation occurs in the bulk. As 
shown in [5], this model enjoys strong chaotic properties. 
Its natural invariant measure is fractal, with one posi- 
tive and one negative Lyapunov exponents, whose sum 
is negative, and identified as minus the entropy produc- 
tion rate [6']. The comparison with the corresponding 
phenomenological expression provides a relation between 
the phase space contraction rate and conductivity. 

The GIK dynamics can actually be given a Hamilto- 



nian formulation [7 
general formalism 



This result is contained in a more 
, by which one identifies the GIK 
trajectories with the geodesies of a torsion free connec- 
tion, called the Weyl connection. In this formalism, one 
shows the GIKLG is equivalent, by a conformal trans- 
formation, to a distorted billiard whose trajectories are 
straight lines, referred to as W-flow. In this geometry, 
the cells are stretched and trajectories accelerated in the 
direction of the external field. A detailed analysis of the 
W-flow Q reveals that the cells' stretching is the essen- 
tial mechanism to driving the non-equilibrium stationary 
state. 

In this letter, we investigate the connection between 




FIG. 1: (Color online) The self-similar Lorentz channel bil- 
liard and a trajectory. 



the chaotic dynamics and non-equilibrium thermody- 
namics of a self-similar billiard chain, such as defined in 
[lo| . This billiard chain can be thought of as a simplified 
version of W-fiow described above; we retain the sretched 
geometry of the cells, but get rid of the bulk dissipation. 

A self-similar billiard is thus a one-dimensional spa- 
tially extended system such as shown in Fig. [1] It consists 
of an infinite collection of two-dimensional cells, identi- 
cal in shapes, but differing in sizes, which are scaled and 
glued together along a horizontal axis. Each cell con- 
sists of convex obstacles upon which point particles col- 
lide elastically, independently of one another. The left 
and right borders of every cell are open and scaled by 
a constant ratio which we denote by jj.. The system is 
so constructed that /x = 1 corresponds to the Lorentz 
channel which exhibits diffusion, but no bias. When 
/i 7^ 1, the system sustains a non-equilibrium stationary 
state, characterized by a steady density current from the 
smaller to the larger scales. 

Strictly speaking, this self-similar billiard channel is 
different from the W-fiow, constructed along the lines of 
0. However the two systems have in common the scal- 
ing property between the left and right borders of ev- 
ery cell. An important difference concerns the dynamics 
between collisions. Here we assume free propagation of 
point particles between collision events, which amounts 
to suppressing the thermostating mechanism of the GIK 
dynamics. Another example of a billiard with a similar 
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scaling mechanism was considered in ll| . Other systems 
which bear similarities with this one are the multibaker 
maps with energy [l^ . Both systems are area- preserving 
and energy conserving, with an external driving force 
(here in the form of a geometric constraint) that induces 
a macroscopic current, without the necessity of a ther- 
mostating mechanism. 

In order to specify the stationary properties of the sys- 
tem, we resort to the special flow of the billiard [13| , thus 
substituting the real time dynamics by a discrete map- 
ping from one collision event to the next. This is a usual 
procedure for the study of ergodic properties of billiards, 
which allows us to replace the volume-preserving dynam- 
ics on the extended system by a phase-space contracting 
mapping on a periodic cell. This feature is also common 
to multibaker maps with energy, which can be reduced 
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to dissipative baker maps such as studied in 
[l^ . Moreover the invariant measure has fractal proper- 
ties much like that of the GIKLG J^. 

Following 0, we identify the phase space contraction 
rate with the entropy production rate, and take the limit 
/i ^ 1 so as to compare this expression to its thermody- 
namic counterpart. We thus infer a relation between the 
current of the near-equilibrium system (/i 7^ 1) and the 
diffusion coefhcient of the equilibrium system {fi = 1), 
much like a fluctuation-dissipation theorem. This result 
is confirmed by our numerical computations. 

Going back to the definition of the model, the reference 
cell, with index i = 0, is represented in Fig. [2ja). It is a 
region defined by the exterior of five disks, four of which 
are centered at the corners of the cell, and one located 
at the center. The dissymmetry between the left- and 
the right-hand sides depends on the scaling parameter fi. 
The other tunable parameter is the ratio R/d, between 
the center disk's radius R and the cell's horizontal width 
d. The radii of the outer disks are fixed to R/ ^JJi (to the 
left) and R^/Ji (to the right). The vertical separations 
between the outer disks is taken to be d\jzjyu to the 
left and dy/3jl to the right, so that the usual Lorentz 
channel is retrieved for /i = 1. We note that restrictions 
must be imposed on the range of permitted values of 
the parameters for the self-similar billiard to share the 
hyperbolicity of the Lorentz channel. The details can be 
found in [10[. We assume these conditions to be met in 
the sequel. 

The whole chain is constructed by adding a cell to the 
right of the reference cell, with index i = 1, identical 
in shape but with all its lengths multiplied by fi, and 
another one to the left, with index i — —1, with all the 
lengths divided by /i. We repeat this construction in 
such a way that in the ith cell, z G Z, all the lengths 
are multiplied by /i*. The resulting billiard chain is so 
constructed that the mirror symmetry with respect to 
the transformation fi ^ 1/fi remains. 

Consider a particle moving inside the billiard with ve- 
locity V. Figure [T] shows such a trajectory. As the parti- 
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FIG. 2: (a) Geometry of the reference cell, (b) Sketch of the 
boundary conditions in the periodic cell, as given by Eqs. ^ 



cle moves from one cell to a neighboring one, the length 
scales change by a factor /i. Equivalently we can rescale 
the velocity by a factor and keep the length scales 
unchanged. For both transformations the characteristic 
time between successive collisions with the walls changes 
by a factor fj,. 

We can therefore analyze the dynamics on the ex- 
tended self-similar billiard in terms of the dynamics in 
a periodic cell by applying proper boundary conditions 
and rescaling the vertical coordinate and velocity, as the 
particle exits either 
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f — 


(d, h) 


V ~ 




r ~ 


iO,h) 


V — 





(o,Vm), 

{Vx/fi,Vy/n). 

{d, hfi), 
{vxH,Vyn). 



(1) 
(2) 



Here h is the y coordinate of the trajectory as it crosses 
the cell's left or right boundaries. An example is depicted 
in Fig. [2tb) . These conditions are analogous to periodic 
boundary conditions for the usual Lorentz channel bil- 
liard, but for the self-similar structure parametrized by 

Billiards are instances of so-called special flows |13| . 
which is to say their time evolutions can be specified 
by the mapping of the coordinates between successive 
collision events, together with the time interval between 
them. This kind of representation is particularly useful 
to the evolution of self-similar billiards, where the ab- 
sence of an overall time scale in the extended system 
poses difficulties in the definition of proper time aver- 
ages. Indeed point particles typically move towards the 
direction of increasing cell sizes, thereby decreasing their 
collision frequency. Such ambiguities become irrelevant 
when the evolution is considered from one collision event 
to the next. 

The particle's coordinates from one collision with the 
walls to the next are specified by the Birkhoff map Q 
= (SjVjtu)^ 1-^ = Here the variable s 
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represents the arc-length along the unit cell boundary 
(including the open sides), v is the modulus of the parti- 
cle's velocity and w is the sine of the angle between the 
outgoing velocity and the normal to the cell's boundary. 
This map, together with Eqs. (PH]), provide the map for 
the specification of a trajectory. 

In order to keep track of the correspondence between 
the dynamics on the periodic cell and that on the ex- 
tended lattice, we introduce a new variable /„, which 
takes integer values and labels the cell where the particle 
is located after n iterations of the map. This variable 
changes according to /„+i = /„ + 0(61)7 where we intro- 
duced the jump function, a(^„), which takes the values 
a(Cn) = 1 (resp. —1) if 4'{£,n) has the spatial coordi- 
nate Sn+i on the right (resp. left) open side of the cell, 
otherwise a(^„) = 0. 

The time it takes a trajectory at a (phase space) point 
^ on the boundary to intersect again with the boundary 
of the billiard depends on the speed v as 
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(3) 



where L(^) is the length of the trajectory between inter- 
sections at ^ and (/)(^) of the trajectory with the boundary 
of the unit cell. A complete characterization of the con- 
tinuous time flow is thus given by the variables (^ r) in 
the unit cell, with < r < r(^), a new variable that 
restores the position between collisions. The equivalence 
between the continuous time flow on the extended lattice 
and the Birkhoff map with the extra variable r is further 
discussed in [lO|. 

A remarkable property of the Birkhoff map defined 
above is that it does not preserve phase-space volumes. 
Indeed, as the particle exits to the right and re-enters 
to the left or vice-versa, both s and v coordinates are 
rescaled by fi^^ (resp. /i). Accordingly, probability mea- 
sures evolve under iteration of the density evolution op- 
erator associated to the Birkhoff map towards a unique 
invariant measure (i. e. a probability measure invariant 
under this operator) with fractal properties, characteris- 
tic of a non-equilibrium stationary state. This measure 
has three Lyapunov exponents, two negative and one pos- 
itive, Ai > > A2 > A3, and such that the phase space 
contraction rate cr = —(Ai -|- A2 -f A3) > 0. 

The second Lyapunov exponent, A2, is specific to the 
self-similar geometry of the billiard ; it is due to the con- 
traction of the V coordinate as the particle moves around 
the cell : 



A2 



lim — — log/i < 0, 
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where A/f, = /„ — /q is the lattice displacement vector 
after n iterations, or winding number. If /x > 1, the par- 
ticle moves preferentially to the right, corresponding to 
A/„/n > and the other way around if /i < 1. The 
symmetric case /i = 1 is the usual Lorentz channel with 



only two non-trivial Lyapuonv exponents, i. e. X2 = 0. 
We therefore see that the stationary value of v under 
the Birkhoff map is trivially w„ — > Uqo = 0, which is to 
say particles ever move towards the larger lattice scales, 
where times between collisions keep increasing and colli- 
sions become seldom. 

A2 accounts for one half of the total phase space con- 
traction rate. The contraction of the vertical coordinate 
in Eqs. ([T][21) accounts for the other half, due to Ai and 
A3, 



Ai + A3 = A2 < 0. 



(5) 



These two Lyapunov exponents are related to the com- 
ponents s and ZD of ^ 

The invariant measure therefore has a product struc- 
ture, d"^TO(^) = d^mi(s, ro)dm2(w). When 7^ 1, we 
have dm2{v) = S{v)dv on the one hand, where S is the 
Dirac delta function, and, on the other hand, since Ai > 
and A3 < 0, mi is a fractal , whose Lyapunov dimension, 
here defined by 1 + Ai/jAsj, is somewhere between 1 and 
2. 

The entropy production per collision can be defined 
dynamically as the phase space contraction rate 0, tr = 
— (Ai -|- A2 -l- A3) — — 2A2, which can be computed in terms 
of the winding number, as in Eq. (|4]). 

The winding number can furthermore easily be con- 
nected to the average drift of particles on the extended 
system. As we showed in , the average of the spatial 
coordinate X on the extended system evolves linearly 
in time, with a slope V which defines the drift velocity, 
{X)t = Vt. The winding number and drift velocity are 
equal up to dimensional factors : 
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where d is the reference cell's width, and the aver- 
age time between collisions measured in that same cell. 
The phase space contraction rate can therefore be ex- 
pressed directly in terms of the average drift velocity : 
cr = log /i. This expression can be transposed to a 
phase space contraction rate per unit time in the limit 
/i — > 1. In this limit, V and log/i both scale like /i — 1, so 
that the lowest order contribution to cr is ©(/z — 1)^, with 
T^=i. Divinding cr by the latter defines the (dynamical) 
entropy production per unit time. 
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As we showed in [lo| , the phenomenology of self-similar 
systems such as this one is characterized on the one hand 
by a constant drift from the small to the large scales and, 
on the other, by a mean square displacement which un- 
dergoes a transition from diffusive to ballistic behavior 
at a crossover time which scales like /i'^/^/(/i— 1). There- 
fore, as fj, 1, self-similar systems can be considered as 
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diffusive systems with constant drift, such as modeled by 
a Smoluchowski type equation for the probabihty density 
W, dtW = VdlW - Vd^W. In the stationary state, the 
phenomenological entropy production rate is given 14 1 
by the drift velocity V squared, a term quadratic in ^ — 1, 
divided by the diffusion coefficient V, here computed at 

^i = l. 

Equating Eq. ((T]) with the expression of the thermo- 
dynamic entropy production rate at 0(fi ~ 1)^, we infer 
the following relationship between the velocity drift and 
diffusion coefficient : 



lim 

M^i /i — 1 



2V 



(8) 



which can be viewed as a definition of the mobility, the 
external field being given by 2(/i— 1), and therefore estab- 
lishes the equivalent of a fluctuation-dissipation theorem. 
This relationship is confirmed by the numerical results 
presented in Fig. [31 which were carried out for different 
choices of R/d and four different values of /z close to 1. 
Equation ([6]) was also verified in the numerical scheme. 
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FIG. 3: (Color online) The difference 2V / d - V / - 1) vs. 
R/d is — 1)^, as expected from Eq. ([SJ. Here shown are 
the values /i = 1.02, 1.05, 1.10, 1.15 from bottom to top. The 
straight lines show 2(/i — 1)^. Here V is computed as the 
slope of {X)t vs. t, where the average is computed out of 
many random initial conditions. 

To conclude, self-similar billards are simple mechani- 
cal models of conduction in spatially extended systems 
with volume-preserving dynamics, whereby a geometric 
constraint induces a steady current. As in the conformal 
transformation of the GIKLG 0] , or the multibaker maps 
with energy the essential ingredient for this behav- 
ior is the rapid increase of phase-space volumes, whose 
rate is here given by ^. 

These models lend themselves to further exploring 
the connections between the chaotic properties of non- 
equilibrium dynamical systems and phenomenological 
thermodynamics. The dynamical properties are conve- 
niently analyzed by the means of the Birkhoff map on 



the periodic cell, which specifies the evolution of phase 
points from one collision to the next, here on a three- 
dimensional space. 

Whereas the dynamics on the extended system pre- 
serves phase-space volumes, the Birkhoff map on the pe- 
riodic cell contracts phase-space volumes, with an invari- 
ant measure whose fractal dimension is between 1 and 2. 
A similar mechanism of contraction of phase-space vol- 
umes at the periodic boundaries occurs with the W-flow. 
The phase-space contraction rate can furthermore eas- 
ily be computed in terms of the map's winding number, 
which bears a simple connection to the average velocity 
drift associated to the macroscopic current. The com- 
parison between phase-space contraction rate and phe- 
nomenological entropy production yields an expression 
of the velocity drift near the equilibrium regime. 
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